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Abstract 

We investigate the possibility to use the neutrinos coming from a future galactic supernova explosion to perform 
neutrino oscillation tomography of the Earth's core. We propose to use existing or planned detectors, resulting 
in an additional payoff. Provided that all of the discussed uncertainties can be reduced as expected, we find that 
the average matter densities of the Earth's inner and outer cores could be measured with a precision competitive 
with geophysics. However, since seismic waves are more sensitive to matter density jumps than average matter 
densities, neutrino physics would give partly complementary information. 
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1. Introduction 

In order to obtain more information about the 
interior of the Earth, neutrino tomography has 
been considered as an alternative method to geo- 
physics. There exist, in principle, two differ- 
ent such techniques, neutrino absorption tomog- 
raphy [UEIEIIIIEIEIEIIEIE! and neutrino oscilla- 
tion tomography |lUUllllT t 2lll3lll4| . Neutrino ab- 
sorption tomography, based on the absorption of 
neutrinos in matter, is in some sense similar to 
X-ray tomography and unfortunately faces sev- 
eral problems including the need of extremely 
high energetic neutrino sources, huge detectors, 
and the prerequisite of many baselines. Neu- 
trino oscillation tomography uses the fact that 
neutrino oscillations are influenced by the pres- 
ence of matter ^JEOEl- Neutrino oscillation 
tomography would, in principle, be possible with 
a single baseline, since interference effects pro- 
vide additional information on the matter density 
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profile. However, it requires quite precise knowl- 
edge about the neutrino oscillation parameters 
and stringent bounds on the contribution of non- 
oscillation physics, such as neutrino decay, CPT 
violation, non-standard interactions, sterile neu- 
trinos, etc. Supernovae as neutrino sources are 
especially interesting, since the neutrinos come in 
large numbers from a short burst, which could be 
used to obtain a snapshot of the Earth's interior. 
In addition, compared to solar neutrinos, their 
energy spectrum has a high-energy tail, which is 
more sensitive to Earth matter effects. The influ- 
ence of Earth matter on supernova neutrinos has, 
for example, been studied in Refs. |18U19ll2U] . 

We assume that technologically feasible de- 
tectors exist, such as Super-Kamiokande, SNO, 
Hyper-Kamiokande, and UNO, which are orig- 
inally built for different purposes, but also ca- 
pable to detect supernova neutrinos. We dis- 
cuss the possibility to use the neutrinos coming 
from a future galactic supernova explosion to de- 
termine with the assumed detectors some of the 
measurable quantities describing the structure of 
the Earth's interior. We especially focus on the 
outer and inner core of the Earth, since they are 
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much harder to access with conventional geophys- 
ical methods than the mantle of the Earth. 



2. Geophysical aspects 

In geophysics, the most promising technique to 
access the Earth's interior is to use seismic wave 
propagation (for a summary, see Refs. |21H22j '). 
Especially, the detection of seismic waves pro- 
duced by earthquakes gives valuable information 
on the seismic wave velocity profile of the Earth 
matter. However, the matter density is not di- 
rectly accessible, but indirectly obtained by as- 
sumptions about the equation of state of the 
considered materials. Since seismic S-waves are 
mainly reflected at the mantle-core boundary, in- 
formation on the Earth's core is much harder to 
obtain than on the Earth's mantle. Therefore, 
we will especially focus on the Earth's core in 
this paper. Since reflection and refraction of P- 
waves at transition boundaries with large matter 
density contrasts are quite easy to observe with 
seismic waves, the mantle-core and outer core- 
inner core boundaries can be located with high 
precision as well as the matter density jumps can 
be measured. For example, the matter density 
jump at the outer core-inner core boundary is of- 
ten given by (0.55±0.05) g/cm 3 23 . This is quite 
a difference compared to neutrino physics, since 
neutrino oscillations in matter are especially sen- 
sitive to the average matter densities (on the scale 
of the neutrino oscillation length) indirectly mea- 
sured by the electron density. Thus, neutrino os- 
cillations are less sensitive to local differences, but 
they involve less unknowns from the equation of 
state and could therefore access the absolute mat- 
ter densities instead of the matter density jumps. 

Several issues regarding the Earth's inner core 
are considered to be interesting from a geophys- 
ical point of view. For different indirect reasons 
the inner core is believed to consist mainly of iron 
and it is therefore often called the iron core. First 
of all, the spectral lines in the sunlight indicate 
that the atmosphere of the Sun partly consists of 
iron as a potential material source for the plan- 
ets. A second hint comes from the magnetic held 
of the Earth. After all, there are no convincing 
alternatives. We will sec later that neutrino to- 



mography could directly verify the average mat- 
ter density of an iron core. Further topics rele- 
vant for the inner core structure are: anisotropy, 
heterogeneity, time-dependence, solidity, and ro- 
tation (for a summary, see Ref. HU). However, 
since neutrino oscillations are to a good approxi- 
mation only sensitive to average matter densities 
at long scales, these issues are much harder to 
access with neutrino oscillation tomography. 



3. Core-collapse supernovae as neutrino 
sources 

Core-collapse supernovae represent the evolu- 
tionary end of massive stars with a mass M > 
8-M©, where M Q is the mass of the Sun. In these 
explosions, about 99 % of the liberated gravita- 
tional binding energy, Eb — 3 • 10 53 erg, is carried 
away by neutrinos in roughly equal amounts of 
energy for all flavors in the hrst 10 seconds after 
the onset of the core collapse j^S] • It is widely be- 
lieved that the time-dependent energy spectrum 
of each neutrino species can be approximated by 
a "pinched" Fermi-Dirac distribution |2lfll27ll28j . 
In this work, we assume that the time-integrated 
spectra can also be well approximated by the 
pinched Fermi-Dirac distributions with an effec- 
tive degeneracy parameter rj, i.e., 
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where £^, ot = J L v dt is the total neutrino energy 
of a certain flavor i/, 
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c x-y + l 



dx, 



and E and T denote the neutrino energy and the 
effective temperature, respectively. The relation 
between the average neutrino energy and the ef- 
fective neutrino temperature is given by 
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In this section, we assume for simplicity that rj 
for all neutrino flavors 1 , which means that 



(E v) =mT= 



7tt 4 



F 2 (0) 180C(3) 



T= kT, 



(3) 



where £(x) is the Riemann z-function [C(3) — 
1.20206] and k ~ 3.15137. Furthermore, the time- 
integrated flux of the neutrinos can be expressed 
by 



N° 



(4) 



where D is the distance to the supernova. 

Due to the different trapping processes, the dif- 
ferent neutrino flavors originate in layers of the 
supernova with different temperatures. The elec- 
tron (anti) neutrino flavor is kept in thermal equi- 
librium by (3 processes up to a certain radius 
usually referred to as the "neutrinosphere" , be- 
yond which the neutrinos stream off freely. How- 
ever, the practical absence of muons and taus 
in the supernova core implies that the other two 
neutrino flavors, here collectively denoted by u x 
{i^ii-, Vti S'/ij fy)j interact primarily by less efficient 
neutral-current processes. Therefore, their spec- 
tra are determined at deeper, i.e., hotter, regions. 
In addition, since the content of neutrons is larger 
than that of protons, v e 's escape from outer re- 
gions than i> e 's. This rough picture leads to the 
following hierarchy: (E Ue ) < (E 9c ) < {E Vx ). 
Here v x refers again to both and v T . Typ- 
ical values of the average energies of the time- 
integrated neutrino spectra obtained in simula- 
tions are (E Ue ) ~ 12 MeV, (E Pe ) ~ 15 MeV, and 
(E Vw ) ~ 24 MeV |5HI5T] . However, recent studies 
with an improved treatment of neutrino trans- 
port, microphysics, the inclusion of the nucleon 
bremsstrahlung, and the energy transfer by re- 
coils, find somewhat smaller differences between 
the v e and v x spectra |32II33| . 

In the following, we assume a future galactic 
supernova explosion at a typical distance of D = 
10 kpc, with a binding energy of Et, = 3- 10 53 



erg 



1 In the analysis, we will also show some results for r\v^ = 3 
and r\o = 1 1291 . The extension of the formulas to these 
cases are quite straightforward, though making them more 
complicated. 
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Table 1 

Our standard scenarios for supernova parameters, 
where in all cases (Ep B ) = 15 MeV. 



and a total energy of Elf = Eff = ££°7£, 
where £ parameterizes a possible deviation from 
energy equipartition We also assume that 

the fluxes of , v T , v^, and v r are identical 
and we fix (Ep c ) to 15 MeV. Due to the cur- 
rent lack of a standard picture of core-collapse 
supernovae we consider for our analysis six sce- 
narios for rj = with different values of the pa- 
rameters t e = (E Pns ) I (E P J € {1.2, 1.4} and 
£ G {0.5, 1, 1.2}, as shown in Tab. n In addition, 
we study four scenarios with rjp c — 3 and rjp^ = 1 
as well as t e € {1.2, 1.4} and £ e {1, 1.2}. 

As far as the neutrino detection is concerned, 
we only analyze the charged-current reaction v e + 
p —> e + + n, since this reaction yields the largest 
number of events (around 8000 in the Super- 
Kamiokande detector in the case of a galactic su- 
pernova). Therefore, we shall concentrate on the 
study of the propagation of antineutrinos from a 
supernova to detectors on the Earth. The cross 
section details for the reaction D e + p — > e + + n 
can be found and are discussed in Ref. [3*5] , 

4. From neutrino production to neutrino 
detection 

In general, neutrino propagation from a source 
to a detector is described by an evolution operator 
on the form 



U = U(L) = e- i So^')^\ 



(5) 
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where the exponential function is time-ordered, 
ffl = J$?(L) is the total Hamiltonian and L is the 
neutrino (traveling) path length, i.e., the base- 
line length. The Hamiltonian is usually given 
either in the flavor basis or in the mass basis. 
In the flavor basis, the total Hamiltonian reads 
Jf f (L) = UHmlJ- 1 + A(£)diag (1,0,0), where 
H m = diag (Ei, E2, E3) is the free Hamiltonian 
in the mass basis, U = (U aa ) is the leptonic mix- 
ing matrix, and 

A = A(L) = ±V2G F —p(L) (6) 
m N 

is the mass density parameter related to the mat- 
ter density p = p(L). Here E a = m 2 a /(2E) (a = 
1,2,3), G F ^ 1.16639 • 10" 23 eV" 2 is the Fermi 
weak coupling constant, Y e is the average number 
of electrons per nucleon 2 , and win — 939 MeV 
is the nucleon mass. The sign depends on the 
presence of neutrinos (+) or antineutrinos (— ). 
Furthermore, m a is the mass of the ath mass 
eigenstate u a and E is the neutrino energy. In 
order to obtain the neutrino oscillation transition 
probabilities, we need to calculate the matrix ele- 
ments of the evolution operator in the flavor basis, 
take the absolute values of these, and then square 
them. The neutrino oscillation probability ampli- 
tude from a neutrino flavor u a to a neutrino flavor 
up is defined as 

A a/3 = (up\Uf(L)\u a ), a,/3 = e,p,r, (7) 

where 11 j is the total evolution operator in the 
flavor basis. Then, the neutrino oscillation tran- 
sition probability for u a — > up is given by 

P a p=\A a p\ 2 , a,(3 = e,p,r. (8) 

The initial neutrino fluxes arise from the cen- 
tral part of the supernova, where the matter den- 
sity is of the order of about 10 12 g/cm 3 . For such 
high matter densities one can infer from the ex- 
pression of the Hamiltonian Jtf'f that the matter 
mass eigenstates, u™ (a = 1,2,3), coincide with 
the flavor states, u a (a — e, p, r), up to a rotation 
between u^ and u r . Thus, in the case of normal 
mass hierarchy, m\ < to 2 <C OT3, one has 

u? = u e , u? = u>^, u? = u' T , (9) 

2 In the Earth: Y e ~ ~. 



where u' and uL are the rotated states. Tliere- 
^i T 

fore, one can assume that the original fluxes of 
the matter mass eigenstates are 3 

3° r = $° e , $° r = ^ , 4>o„ = $« x . (10) 

Since we are assuming that u^ and u T have the 
same fluxes, the neutrino transitions are deter- 
mined by the mixings of the u e only, i.e., by 
U ea |37| . Moreover, under the assumption of 
normal mass hierarchy, antineutrinos do not un- 
dergo any resonant conversion, which means that 
the small mixing angle #13 is suppressed in mat- 
ter and the u e <-> D3 transitions are negligible. 
One consequence is that u™ propagates adia- 
batically and leaves the supernova as u 3 . The 
propagation of the other two states depends on 
the parameters of the solution to the solar neu- 
trino problem and it may be adiabatic or non- 
adiabatic [SHUSH]- In particular, we will focus 
on the Mikheyev-Smirnov-Wolfenstein (MSW) 
[1511161 IT?] large mixing angle (LMA) solution, 
since it is by far the most favored one. For the pa- 
rameters within such a region the neutrino evolu- 
tion through the supernova (SN) envelope is adia- 
batic. Thus, u e will leave the supernova as U\, u' 
as U2, and u' T as u 3 . Finally, the measured fluxes 
of supernova neutrinos at a detector on the Earth 
are 

3 

®u a = 5I PaQ$ °a' a = e,p,r, (11) 

a=l 

where $° (a — 1,2,3) are the initial supernova 
neutrino fluxes. 

For neutrino propagation from a supernova (SN 
core) to a detector at the Earth (see Fig. ^) we 
have the probability amplitudes 

A a0 = (vp\Uf\L)\v a ) 

= (v \U®{L®)U™{L v ^)U w {L S k)K), 
(12) 

3 Any rotation between 9^ and P T does not affect the cor- 
responding total mass eigenstate contents, because they 
have the same fluxes, as discussed in the last section. For 
an analysis taking into account possible differences in the 
fluxes, see Ref. 1361 . In addition, an argument against neu- 
trino oscillations between and v T on their way to the 
Earth will be given at the end of this section. 
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Earth 



SN envelope 



SN core 




Figure 1. The propagation of neutrinos from a supernova (SN) to the Earth. The detector D\ can be 
anywhere on the Earth's side towards the SN, whereas the detector D 2 should be in the shadow of the 
Earth's core. 



where U SN , U vac , and U® are the evolution oper- 
ators in the supernova (from SN core to SN enve- 
lope), in vacuum, and in the Earth, respectively, 
and Lsn, ivac, and L§ are the corresponding 
baseline lengths. Note that the operators U SN , 
U vac , and U® in general do not commute. Us- 
ing the completeness relation, one can write the 
probability amplitudes as 



A aP = ^(^|W®(L®)Z^ ac (L vac )K) 



u a \U^(L SN )\u a ) 



(13) 



Since we have seen that for adiabatic transitions 
the supernova neutrinos leave the supernova (SN 
envelope) as neutrino mass eigenstates v ai i.e., 
Va) — 5 aai we can re-define the 
probability amplitudes 



(u a \U®(L E )U™ c {L vac )\v a ) 



(14) 



where the first index is a mass eigenstate index 
(a — 1, 2, 3) and the second index is a flavor state 
index (a — e, /z, r). These "mixed" probability 
amplitudes will completely determine the evolu- 
tion of the neutrinos from a supernova (SN enve- 
lope) to the Earth. Now, there are, in principle, 
two cases for a supernova neutrino to be detected 
at the Earth (see again Fig. [TJ: 



1. The supernova neutrino arrives at the de- 
tector from above, i.e., it does not go 
through the Earth at all (detector Di). 



2. The supernova neutrino goes through the 
Earth's and then arrives at the detector (de- 
tector D 2 ). 



Let us start with the first case. The probability 
amplitude for an initial neutrino mass eigenstate 
v a , where a = 1, 2, 3, to leave the supernova and 
to convert into a flavor state v ai where a = e,fi, r, 
is at the detector D\ formally defined as 



.4 



;^|zr ac (L vac )K). 



(15) 



Note that we assumed that a neutrino mass eigen- 
state v a left the supernova, and therefore, no evo- 
lution operator 14 SN should appear in the above 
equation. Furthermore, since the supernova neu- 
trino does not go through the Earth, there ap- 
pears also no evolution operator U® in this equa- 
tion. Next, since the evolution operator in vac- 
uum U vac is diagonal in the mass basis, we find 
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that Eq. ijT5)l reduces to 4 



= (VaWa) = > Uab(Vb\u a ) 



3 

E< 

6=1 



3 

E 

b=l 



(16) 



i.e., the probability amplitudes are just the ma- 
trix elements of the leptonic mixing matrix. 
Thus, we have = \Aj& \ 2 = \U aa \ 2 , and in- 
serting this into Eq. , we obtain the supernova 
neutrino fluxes at the detector D\ as 



= 1 



e,/i,T. 



(17) 



Assuming as in Eq. IjlOl) that the initial fluxes of 
the second and third mass eigenstates are equal, 
i.e., <I>? 2 = $? 3 , the electron antineutrino flux at 
the detector Di can be written as 



= $2 [\U el \ 2 + f R (|C/ e2 | 2 + \U e3 \ 2 )] . (18) 

Here f R = $ p J$ Pl = is the so-called 

flux ratio, which is plotted for several values of £ 
and te (introduced in the last section) in Fig. |21 
Furthermore, the flux ratio /jj depends on the 
supernova parameters {E Pc ), £, and te only and 
reads for 77 = 5 



fa = 



: V^J + 1 



'EqTe^J'e _|_ 1 



(19) 



where again k ~ 3.15137. Reinserting the proba- 
bilities P^ 1 instead of the probability amplitudes 
in Eq. (|18fl and using the unitarity condition 



p ie + p ?e x + Pi 1 = 1, we find that 



2e 



(20) 



which means that the flux of electron antineutri- 
nos at the detector D\ is only depending on the 

4 The neutrino flavor states are defined as follows: \v a ) = 
Eo=l UaaWa) ( a = e ,^! T )> which implies that (u a \ = 
5^^ =1 Uaa{va\. Here the Uaa's are the matrix elements of 
the leptonic mixing matrix U. 

5 For r\ ^ this equation would be slightly more compli- 
cated, but could be easily obtained from Eq. IT1 . 
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Figure 2. The flux ratio fa as a function of the 
neutrino energy E for different values of the pa- 
rameters £ and te and for rj — 0, partly corre- 
sponding to the scenarios in Tab. ^ 



initial neutrino flux , the transition probabil- 
ity P-^ 1 = |Z7ei| 2 , and the nux ratio Jr. 

Next, let us discuss the second case. Again us- 
ing the fact that the evolution operator in vacuum 
^vac j g d} a g 0na j m ^ ne ma ss basis, we find that 

A% = <^®(L e )zr ac (L vac )k> 
3 

= ("« ( L © ) I ^) (** I ( L vac ) I ) 

6=1 
3 

= ^(z/ Q |W e (L e )K)<U 

6=1 

= (^W®(L ffi )K>=^. (21) 

Similar to the first case, we obtain the supernova 
neutrino fluxes at the detector D2 as 



< 2 = E 

0=1 



U® | 2 $° 



a=e,n, t, 



(22) 



which for the electron antineutrino flux at the 
detector D2 can be written as 



D 2 



W 1 + /^ 



(23) 



This means that the flux of electron antineutri- 
nos at the detector D 2 depends only on the ini- 
tial neutrino flux , the transition probability 

p ie = \ A ?e\ 2 > and thc nux ratio fR- 
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Now, we want to determine the neutrino oscil- 
lation transition probabilities. Using the evolu- 
tion operator method developed in Ref. 0U|, the 
evolution operator in the Earth, which we will as- 
sume to consist of N different (constant) matter 
density layers, is given by 



x U®{L 2] A 2 )U®{L X ;A X ) 



N 



(24) 



fc=i 



where U®(L k ] A k ) = e - ijr ( Afc ) Lfc is the evolution 
operator in the fcth layer with constant matter 
density and = Ylk=i ^k- 6 Here L k and A k 
are the thickness and matter density parameters 
of the fcth matter density layer, respectively. Note 
again that the evolution operators in the differ- 
ent layers normally do not commute. Inserting 
Eq. l|2"4"|l into Eq. (j2"T|) and the result thereof into 
p aa = \ A ai\ 2 , we finally obtain 



pD 2 _ 



N 



Y[U®(L k ;A k )\v a ) 



k=l 



(25) 



which is our main formula for the neutrino oscil- 
lation transition probabilities from a supernova 
to the detector D 2 . Thus, inserting a — 1 and 
a = e into Eq. 1251) , we find the probability P^ 2 , 
which can then be inserted into Eq. I|23|) . 

In Fig. |2| we show for the scenario SI in Tab.^ 
the different binned energy spectra at Super- 
Kamiokande-like detectors, where the solid curve 
represents the detector D 2 and the dashed curve 
the detector D x . One can easily see that mat- 
ter effects are largest for energies above about 
25 MeV. For a detailed discussion of Earth matter 
effects of supernova neutrinos, see Refs. [18119120] 
and references therein. In summary, for our cho- 
sen values of the neutrino oscillation parameters 
(MSW LMA solution, maximal atmospheric mix- 
ing, and normal mass hierarchy) the relative size 
of the Earth matter effects increases with energy 



"Similar applications of the evolution operator for prop- 
agation of neutrinos in matter consisting of two density 
layers using two neutrino flavors have been discussed in 
Refs. I41I42I43I . 
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Figure 3. The charged-current events per en- 
ergy bin from electron antineutrinos in an Super- 
Kamiokande-like detector for a realistic energy 
resolution 01]. The solid curve shows the spec- 
trum (scenario SI in Tab.^l including the Earth 
matter effects for a baseline of about 12,700 km 
through the Earth's core, whereas the dashed 
curve shows the spectrum without Earth mat- 
ter effects. The differences appear to be small, 
but are nevertheless statistically significant, i.e., 
Ax 2 — 35. Assuming that the matter effects can 
(approximately) be described by two parameters, 
such as the average mantle and core densities, 
this corresponds to at least a 5er effect. To estab- 
lish the Earth matter effects at a 3a confidence 
level, one would need at least 11 kt fiducial mass 
detectors. For the neutrino oscillation parame- 
ters we used 12 = 32.9°, 6 13 = 5°, 9 23 = 45°, 
Amfi = 5.0 • 10" 5 eV 2 , Am| 2 
and Sqp — 0, i.e., a normal mass hierarchy. 



2.5 ■ 10" 3 eV 2 , 



and can be seen as oscillatory modulation of the 
energy spectrum. For small energies, however, 
this modulation oscillates too fast to be resolved, 
i.e., it averages out. For large energies about 
25 MeV, the frequency becomes smaller and this 
modulation can be resolved. In Fig. we used 
the scenario SI from Tab. ^ making it possible 
that the fluxes of v x dominate above the critical 
energy around 25 MeV. This results in a nega- 
tive modulation of the electron antineutrino spec- 
trum, i.e., regeneration effects of and v T , as it 
is shown in Ref. |18| . However, a general sup- 
pression of the fluxes of the v x is possible for a 
value £ < 1, which means that the modulation 
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can be positive. From Fig. [3 it is also interesting 
to observe that solar neutrinos show much smaller 
Earth matter effects (day- night), since the spec- 
trum is cut off far below 25MeV. Thus, espe- 
cially the high-energy tail in the supernova neu- 
trino spectrum can make this application possible 
compared with solar neutrinos. 

Let us now go back to Fig. which shows 
the flux ratio /r of the v x and v\ fluxes at the 
surface of the Earth for some of the supernova 
parameter scenarios in Tab. With Eqs. QJ, 
()19(l . and l)23|l as well as this figure, we can esti- 
mate the sensitivity for different sets of supernova 
parameters. This equation describes the flux at 
the detector Di and depends on three different 
quantities 3?^, f Rl and P® 2 . First, the flux 4?^ 
can be indirectly determined by the detector Di , 
since it also appears in Eq. 1)21) [) as the envelop- 
ing function, i.e., both fluxes "I"? 1 and <I>£ 2 are 
directly proportional to $° . Second, in general, 
the flux ratio depends on the supernova pa- 
rameters £ and te as well as the 77's, which can, 
up to a certain precision, also be reconstructed 
from the spectrum at the detector D\ (c/., dis- 
cussion in Sec. [7| and Refs. 45.:46J). Note that 
it could be directly measured if one were also 
able to detect flavors other than the electron an- 
tineutrino, and the supernova parameters would 
completely drop out (c/., Eq. 1)17(1 ). Third, the 
transition probability P® 2 contains the informa- 
tion about the Earth matter and is usually quite 
large. Thus, the ratio (1 - P° 2 )/P° 2 in Eq. (23 
is rather sensitive to changes in the Earth matter 
effects. Since this factor is multiplied with fn, the 
energy-dependent flux ratio can enhance or sup- 
press it. Finally, we have also noticed above that 
the (relative) Earth matter effects are increasing 
with energy. For the supernova parameters we 
can then distinguish four different cases, where 
some of those can also be found in Fig. [21 

1. £ = 1, T E = 1, r\u e = = (energy 
equipartition and equal temperatures for all 
flavors): The flux ratio /jj is equal to unity 
(c/., Eq. H19JI ). Then, the neutrino transi- 
tion probabilities in Eq. Q23jl drop out and 
we cannot use the supernova neutrinos for 
Earth matter effects. 



2. £ = 1, t e > 1, ?7p e = »?p m = (energy 
equipartition and a lower temperature for 
v e than for v x ): The flux ratio fn is en- 
hanced for large energies, where Earth mat- 
ter effects are large. The larger te is, the 
larger becomes this effect. Thus, the sce- 
nario S2 in Tab. ^ performs worse than the 
scenario SI. 

3. ( < 1 (more electron antineutrinos pro- 
duced than the other two flavors): The 
flux ratio fu oc £ is suppressed in general. 
Therefore, the scenarios S3 and S4 are not 
as good for our application as the scenario 
SI. 

4. te > 1 and £ > 1 or ?7 Po ^ 77^ 
(more muon/tau antineutrinos produced 
than electron antineutrinos or different de- 
generacy parameters for the different fla- 
vors): The flux ratio fu becomes even 
steeper than the one for SI, reflecting the 
more different behavior of the different fla- 
vors. Thus, we expect that the scenarios S5 
to S10 in Tab. Hperfo™ 1 much better than 
the scenario SI, while the optimal scenario 
should be S9 having the largest £ and te- 

Thus, the larger £ and te are and the more differ- 
ent r)p e and r)p are, the better our proposed ap- 
plication should work. However, the actual values 
of these parameters will not be known for sure be- 
fore the next supernova explodes. Therefore, we 
choose further on the scenario SI for reference, 
which has quite plausible parameter values and 
will turn out to be a possible model that per- 
fectly illustrates our procedure. In addition, we 
will summarize the performance of the other sce- 
narios in short form. 

For our application we assume that neutrino 
mass eigenstates arrive at the Earth and no neu- 
trino oscillations take place between the super- 
nova envelope and the surface of the Earth. This 
can either be justified by the adiabacity condition 
for the propagation within the supernova, making 
mass eigenstates emerge from it, or by decoher- 
ence of neutrino oscillations between the surface 
of the supernova and the Earth. In both cases, 
Eq. (|13fl can be split into two independent factors 
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without interference terms. The issue of wave 
packet decoherence has, for example, been ad- 
dressed in Refs. |T71I%1115T15T1I52] . It has been 
found that neutrino oscillations vanish for neu- 
trino propagation over distances much larger than 
the coherence length of the neutrino oscillations. 
This means that for L > L c a f l oc aE 2 /Am 2 ab the 
L/_E-dependent interference terms produced by 
the superposition of the mass eigenstates m a and 
uib in the neutrino oscillation formulas are aver- 
aged out. The quantity a corresponds to a wave 
packet width determined by the production and 
detection processes 48 , 50 . Since for supernova 
neutrinos the distance of the propagation is espe- 
cially large, it is plausible to assume that this 
averaging takes place and neutrino oscillations 
vanish by natural decoherence. In other words, 
the different group velocities of the wave packets 
of different mass eigenstates combine with dis- 
persive effects such that the overlap of the mass 
eigenstates is gradually reduced to zero by a fac- 
tor of exp (— [l/L c °^] 2 ) in the neutrino oscillation 
formulas. Hence, for L > L c °^ the mass eigen- 
states arrive separately and the coherent transi- 
tion amplitude in Eq. I|12|) is split up into two 
parts to be summed over incoherently (see, e.g., 
Ref. 51 ), having the same effect as the adia- 
batic transition within the supernova separating 
the flavor states into mass eigenstates. Thus, it is 
reasonable to assume that mass eigenstates arrive 
at the surface of the Earth even for non-adiabatic 
transitions within the supernova. 

Finally, in our numerical analysis, we assume 
Supcr-Kamiokandc-like water-Cherenkov detec- 
tors, i.e., a 32 kt fiducial mass (for supernova 
neutrinos) Super-Kamiokande detector and a 
1 Mt fiducial mass (for solar neutrinos) Hyper- 
Kamiokande detector. We choose 30 energy bins 
between the threshold energy 5 MeV and 70 MeV, 
since above 70 MeV the number of events is rather 
low. Furthermore, we assume a realistic en- 
ergy resolution implemented with Gaussian en- 
ergy smearing with a smearing width propor- 
tional to \f~E such that we have an energy res- 
olution of about 15 % at 10 MeV 



5. A neutrino oscillation tomography 
model 

We now introduce and discuss a simple model 
used for supernova neutrino tomography. As 
shown in Fig. we assume at least two base- 
lines ending at detectors with similar statistics 
and systematics, such as Super-Kamiokande-like 
water-Cherenkov detectors. In order to measure 



' SN neutrinos 
(mass eigenstates) 



Neutrinos 
start to 
oscillate 



Mantle 
Core 

D 2 \ /\ 

\ / Inner core 



Figure 4. The minimal required setup for super- 
nova neutrino tomography with at least two base- 
lines, one of which is ending at the surface of the 
Earth at the detector Di, the other one is go- 
ing through the Earth's core (or even the Earth's 
inner core) to the detector D 2 . In general, we 
assume that the neutrinos arrive as mass eigen- 
states at the detector D\ and start to oscillate 
when they enter the Earth's interior. 



the reference spectrum of the supernova neutri- 
nos, the neutrinos detected at the detector D\ 
must not cross the Earth. If we want to obtain 
information on the Earth's core, then the second 
detector D 2 needs to observe the supernova neu- 
trinos with a baseline crossing the Earth's core 
with a sufficient length. The probability to have 
such a configuration depends upon the location of 
the supernova in our Galaxy, the time of the day 
at which the burst arrives at the Earth, and the 
position of the detectors D\ and D 2 themselves. 
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For example, if we consider a supernova located 
in the galactic center and the detectors placed in 
Japan and the USA 7 , then, following Ref. |18| . 
one obtains that the probability to find the re- 
quired setup is around 25 % and 10 % along a 
day for neutrinos crossing the core and the inner 
core, respectively. Thus, although this configura- 
tion would not work out for many locations of a 
supernova in the galactic plane, it would be quite 
likely for a supernova being in the region of the 
galactic center. 

Further on, we assume that the detector D\ is 
at least as good as the detector D2, which means 
that the statistics is limited by D2, and D\ mea- 
sures the reference flux <frp* with sufficient pre- 
cision. For the neutrino oscillation parameters 
we choose 6 12 = 32.9°, 6 13 = 5°, (9 23 = 45°, 
Am 2 ., = 5.0 • 10~ 5 eV 2 , Am| 2 = 2.5 • 10" 3 eV 2 , 

a normal mass 



'21 — 

and (Sep 

hierarchy. Matter effects on supernova neutrinos 
in the Earth are discussed in detail in Ref. |18|. 
where it is also demonstrated that Earth matter 
effects would be suppressed for solutions other 
than the MSW LMA solution. In addition, for 
detecting antineutrinos the Earth matter effects 
would only be large enough for our application 
with an inverted mass hierarchy if |J7 e 3| 2 is larger 
than about 10~ 5 jTHj. Thus, we specialize on the 
(not unlikely) case of the MSW LMA solution 
with a normal mass hierarchy in order to be able 
to observe Earth matter effects with a sufficient 
precision. For our application the dominant neu- 
trino oscillation parameters are the solar neutrino 
oscillation parameters Am 2 ., and O12, as well as 
the matter effects depend on sin 2 29 Later, we 
will estimate the precision with which we need 
to know these parameters and we will test the 
influence of the size of sin 2 26*13. Note that we as- 
sume mass eigenstates arriving at the surface of 
the Earth. Therefore, if D\ and D2 were identi- 
cal detectors, a direct comparison of their energy 
spectra would verify the existence of Earth mat- 
ter effects immediately. For detectors of different 
types fits of the energy spectra would supply sim- 
ilar information in an indirect way. 



For the modeling of the matter density pro- 
file we choose the Preliminary Reference Earth 
Model (PREM) profile (531201, as it is shown 
in Fig. [SI and approximate it by layers of con- 
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7 For instance, one possible site for the UNO proposal 
would be the WIPP in Carlsbad, New Mexico, USA 



Figure 5. The model for the Earth's matter den- 
sity profile used in the calculations (step function) 
and the PREM matter density profile as function 
of the path length along the baseline shown in 
Fig. 2] The quantities which we are interested in 
are the average outer core matter density p and 
the matter density jump Ap between average in- 
ner and outer core densities, as it is shown in this 
figure. 



stant average matter densities. A baseline with a 
maximum length of twice the Earth radius then 
crosses the following layers in this order: mantle, 
outer core, inner core, outer core, mantle. Since 
substantial knowledge is provided by geophysics 
about the Earth's mantle, we assume its matter 
density to be known with a sufficient precision. 
For such a baseline the interesting quantities to 
measure are the average outer core matter den- 
sity p and the matter density jump Ap to the 
average inner core matter density, as illustrated 
in Fig. This is slightly different to what is 
known from seismic wave geophysics, since there 
the density jumps of the actual matter densities 
at the mantle-core and outer-inner core bound- 
aries are better known. However, since neutrino 
oscillations are not sensitive to matter densities 
at individual points, but essentially to the inte- 
gral of the matter density and the length scale 
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the neutrinos are traveling through |6T] , they are 
more appropriate to measure average matter den- 
sities. 

The introduced model allows to estimate what 
could be learned from the neutrinos of a super- 
nova explosion about the Earth's interior. The 
actual situation, however, such as the number of 
detectors with baselines through the Earth, their 
size, the knowledge on the neutrino oscillation pa- 
rameters, systematics, etc., can only be imple- 
mented after the next observed supernova explo- 
sion. Our discussion here serves only the purpose 
of demonstrating that such studies are feasible. 

6. Results 

Based on the modeling in the last section, we 
present results, which could be provided by a sin- 
gle supernova. Our analysis is performed with a 
standard \ 2 technique using [JJ21 

x ^2Y j (xf-x l + x l \og^-X ( 26 ) 

i=l ^ X i ' 

where n is the number of energy bins, a:£ e£ is the 
reference event rate in the ith bin for the true 
parameters, and Xi is the measured/varied event 
rate in the ith bin. The errors quoted are read 
off at the 2a confidence level - depending on the 
problem for one or two degrees of freedom. For 
two degrees of freedom we also take into account 
the two-parameter correlations. However, we as- 
sume that the effects of the systematics are neg- 
ligible, i.e., the systematical errors are not larger 
than the statistical errors and the systematics is 
well understood. This assumption should be rea- 
sonable at the time when this application could 
become relevant. 

The most likely case to observe the Earth's core 
with supernova neutrinos are baselines crossing 
the Earth's outer core, but not the Earth's in- 
ner core, i.e., baselines between about 10,670km 
and 12,510 km. Assuming the mantle proper- 
ties to be known quite well from geophysics, one 
may then measure the average (outer) core mat- 
ter density (about 11.4 g/cm 3 ). As a result of the 
analysis for the scenario SI in Tab. ^ h turns 
out that one could measure this core matter den- 
sity with a baseline just touching the inner core 



(L ~ 12,510 km) with about 9 % precision with a 
Super-Kamiokande-like detector and 1.3 % preci- 
sion with a Hyper-Kamiokande-like detector. We 
find a rather small dependence on the baseline 
as long as L > 11,250 km. For L ~ 11,250 km we 
still find precisions of 16.5 % (Super-Kamiokande) 
and 2.5 % (Hyper-Kamiokande), which are grow- 
ing with an increasing baseline length. The rea- 
son for the weak baseline dependence comes from 
the geometry of the problem, which demonstrates 
that above a certain threshold the traversed core 
fraction is rather large, but hardly changes if the 
total baseline length is further increased. Thus, 
we will further on not discuss the baseline depen- 
dence, since it turns out to be negligible close to 
the maximal traversed core/inner core distance. 

A somewhat more sophisticated application is 
the combined measurement of the outer and in- 
ner core matter densities in the two-parameter 
model introduced in the last section, i.e., Figs. 01 
and 03 In Fig. the results of this analysis 
without systematical errors and uncertainties are 
shown in the p-Ap-plane for Super-Kamiokande- 
and Hyper-Kamiokande-like detectors for the sce- 
nario SI in Tab.n One can read off the (absolute) 
errors at the 2a confidence level such as shown 
in this figure. A Super-Kamiokande-like detector 
cannot verify the inner core, which can be seen 
in the contours crossing the inner core sensitivity 
line, i.e., the line Ap = 0. Furthermore, the 3a 
contour is rather extensive and cut-off in the up- 
per left corner. Therefore, it is not well-suited for 
density measurements of the inner core. However, 
a Hyper-Kamiokande-like detector can clearly ob- 
serve and verify the inner core even at the 99 % 
confidence level. Furthermore, quite precise mea- 
surements of p and Ap are possible. The relative 
error for p is about 3.1 % and for Ap about 59 % 
(2a confidence level). 

As we have indicated in Sec. 01 the supernova 
parameter scenarios in Tab. pother than the sce- 
nario SI perform somewhat worse in the mea- 
surement of the parameters p and Ap. In addi- 
tion, degenerate solutions appear at the 2a confi- 
dence level, i.e., different solutions in the param- 
eter space can be fit to the results of the mea- 
surement at the considered confidence level. In 
Tab. [21 we show the errors for the measurement 



12 




Figure 6. The la, 2a, and 3cr contours of the x 2 -function for a measurement of p and Ap for a Super- 
Kamiokande-like (left-hand plot) and Hyper-Kamiokande-like (right-hand plot) detector and the scenario 
SI in Tab. H without uncertainties and systematical errors. For the detectors, we use the ones described 
at the end of Sec.0] The errors from statistics and correlations are read off at the arrows in the figures. 
In order to find evidence for the existence of the inner core (the iron core), the contour of the respective 
confidence level must not cut the inner core sensitivity line corresponding to Ap = 0. This line can only 
be seen in the left-hand plot. In addition, in the left-hand plot, degenerate solutions are, in principle, 
present at the 2a confidence level, which are, however, out of the plot range. 



of p and Ap for the Super-Kamiokande detector. 
It clearly reflects the behavior expected in Sec. 0] 
and it demonstrates that large values of £ and te 
as well as different values for the 77's of the differ- 
ent flavors give the best results. It also shows that 
the scenario SI is one of the the most conserva- 
tive choices of the scenarios for which the tomog- 
raphy application would work. In addition, note 
that degenerate solutions are in many solutions 
present due to the energy resolution function of 
the detector. It turns out that a somewhat higher 
(Ep c ) can improve the results, because it supports 
the high-energy tail in the spectrum where matter 
effects are largest. 

The results from this measurement cannot be 
directly compared with the geophysical results, 
because in seismic wave geophysics the matter 
density jumps are much easier accessible than 
the average matter densities. For example, the 
matter density jump between the outer and inner 



cores is believed to be about (0.55 ± 0.05) g/cm 3 
(see, e.g., Refs. |55IISll65ll66p . Translated to the 
2a confidence level, this corresponds to the same 
order of magnitude as the Hyper-Kamiokande 
measurement of about 59 % precision. However, 
the average matter density is much harder to ac- 
cess in geophysics and can only be estimated by 
the long periodic seismic eigenmodes with uncer- 
tainties increasing with depth |2.'ij . The precision 
on the average matter density of about 3 % from 
neutrino physics as well as the measurement of 
the average matter density jump Ap could help 
to understand and complement the geophysical 
information. 

Finally, one could consider more than one base- 
line. If more than one detector observes a super- 
nova through the Earth's core or the inner core, 
then the potential of this technique would be im- 
proved. However, the probability for a single de- 
tector to have a baseline through the core is al- 
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Scenario 










Sp 


S(Ap) 


Degs 


SI 


1.0 


1.4 








2.7 


8.2 


Yes 


S2 


1.0 


1.2 








6.7 


17.5 


Yes 


S3 


0.5 


1.4 








> 8 


> 21 


Yes 


S4 


0.5 


1.2 








> 12 


£ 30 


Yes 


S5 


1.2 


1.4 








2.3 


7.0 


No 


SG 


1.2 


1.2 








5.0 


12.5 


Yes 


S7 


1.0 


1.4 


3 


1 


2.0 


5.9 


No 


S8 


1.0 


1.2 


3 


1 


3.7 


10.3 


Yes 


S9 


1.2 


1.4 


3 


1 


1.7 


5.3 


No 


S10 


1.2 


1.2 


3 


1 


3.1 


9.1 


Yes 



Table 2 

The different supernova parameter scenarios from Tab.^ an d the absolute errors 5p and S(Ap) (in g/cm 3 ) 
for the measurement of p and Ap, respectively, at the 2a confidence level with a Super-Kamiokande-like 
detector. In addition, the appearance of degenerate solutions (Degs) at the 2a confidence level is indicated 
(for Super-Kamiokande) . 



ready quite low, which means that more detec- 
tors would mainly increase the probability that 
one has an appropriate baseline. Thus, we have 
in this paper focused on the case of one baseline 
through the Earth. If really more than one base- 
line went through the core, then the statistics of 
the overall measurement would be improved and 
the result could be estimated by a scaling of the 
detector. Having one large detector and one base- 
line corresponds in this application to different 
detectors at similar positions with their fiducial 
masses adding up to the one of the single detec- 
tor. 

7. Uncertainties 

So far, we have taken into account for our anal- 
ysis only the statistical errors and the larger two- 
paramcter-correlation. Once illustrated the ad- 
vantages of our procedure, we now discuss the 
influence of uncertainties on the measurements 
of the parameters. We do this separately, be- 
cause the sensitivity to the parameter extraction, 
in particular, to the astrophysical parameters, de- 
pends significantly upon the model considered. 

In order to discuss the influence of uncertain- 
ties on the measurements, we estimate the pre- 
cision to which the leading neutrino oscillation 
parameters Am^ and #12 have to be known for 



the measurement and we vary them until we ob- 
serve an effect which is as large as the error of 
the measurement of p or Ap. It turns out that 
these leading parameters have to be known with 
about 1 % precision for the Super-Kamiokande- 
like measurements and with about 0.2 % pre- 
cision for the Hyper-Kamiokande-like measure- 
ments. These precisions should be obtainable on 
the typical timescales of galactic supernova ex- 
plosions. 

The parameter sin 2 2#i3 also has some influence 
on the matter effects. For small values, however, 
the (three-flavor) neutrino oscillations reduce to 
two two-flavor neutrino oscillation schemes one 
describing the solar neutrino oscillations and the 
other one the atmospheric neutrino oscillations. 
We tested the influence of this parameter on our 
applications and we found that it can be safely ne- 
glected as long as sin 2 2613 is not too large. Only 
at the CHOOZ bound minor corrections much 
smaller than the error of our measurements have 
to be performed. However, this bound will be 
reduced in the short term future by planned su- 
perbeam and neutrino factory experiments (for a 
summary of expected boundaries see, for exam- 
ple, Ref. gZi). 

One of the most important uncertainties in this 
measurement comes from measuring the electron 
antineutrinos only. It can be easily seen from 
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Par am. 


SuperK 


HyperK 


Effects 


Z£tot 


-5 % 


~ 1 % - 2 % 


< 1 % 


Z^tot 


- 100 % 


~ 10 % 


large 


Ok.) 


~5 % 


- 1 % - 2 % 


< 5 % - 10 % 




- 10 % 


- 1 % - 2 % 


< 5 % 



Table 3 

The uncertainties on the supernova parameters extracted from the Super-Kamiokande (SuperK) or Hyper- 
Kamiokande (HyperK) measurements from Refs. |45II46| and their estimated effects on our parameter 
measurements as percentage corrections. 



Eq. I|2U|) that the extraction of the fluxes of the 
mass eigenstates from the flux of the electron an- 
tineutrino flavor involves assumptions about the 
supernova parameters, which are entering by the 
flux ratio /r. Thus, if only the electron antineu- 
trinos from the supernova can be measured, then 
the tomography problem will be closely connected 
to the determination of the supernova parame- 
ters at a detector on the surface of the Earth. 
Estimates for the precision of the supernova pa- 
rameter determinations depend on the initial val- 
ues considered. Thus, the more different the v e 
and spectra are formed in a supernova explo- 
sion, the stronger the oscillation effects are, and 
therefore, the better one can measure the neu- 
trino temperatures and luminosities. As an ex- 
ample, we show in Tab. the sensitivities ob- 
tained in Refs. gSIEHj. 8 In addition, this table 
shows the estimated influences on the determina- 
tion of our parameters from a variation of the su- 
pernova parameters in the numerical evaluation. 
The reason for the similar percentage effects for 
Super-Kamiokande and Hyper-Kamiokande is the 
parallel scaling of both problems. It is interesting 
to observe that none of the supernova parameters 
has a strong influence on the tomography prob- 
lem except from the overall energy of the muon 
antineutrinos ££ ot . One can show that this pa- 
rameter has to be known up to about 20 % for 
the Super-Kamiokande measurement and 3 % - 
4 % for the Hyper-Kamiokande measurement in 



8 The results in this reference were obtained for r)p c = 
2.6, r)n^ = 0, and £ = 0.5. In this paper, however, we 
have used slightly different parameter values for which the 
precisions on the supernova parameters in Tab. 131 would 
become somewhat worse by about a factor of two. 



order not to have strong effects on the tomogra- 
phy problem. However, this precision cannot be 
achieved by measuring the electron flavor only. 
Altogether, either needs to be measured by 
different experiments or the fluxes of the muon 
and tau antineutrinos need to be determined si- 
multaneously with the electron antineutrino flux, 
making the supernova parameters entirely drop 
out. This can be seen in Eq. Q17f) . which allows 
the reconstruction of all mass fluxes at the detec- 
tor D\ if all flavor fluxes are measured. 

Another issue in the discussion of uncertainties 
is the parameter Y e in Eq. 10 relating the num- 
ber of electrons to the number of nucleons. Since 
the Earth matter effects in neutrino oscillations 
actually depend on the electron density and not 
on the matter density directly, additional uncer- 
tainties enter in the conversion of these two quan- 
tities by the parameter Y e . We assumed Y e = 0.5 
in our calculations, but for different materials this 
parameter can differ somewhat from this value - 
especially in the inner core. In order to find out 
the material in each matter density layer, one may 
prefer to measure the electron density instead of 
the matter density. However, since in each layer 
these quantities are proportional to each other, 
the problem does not change by using the mat- 
ter density and the electron density can be easily 
calculated. 

So far, we have assumed the average mantle 
matter density to be known exactly. In fact, mat- 
ter density uncertainties up to 5 % have been 
documented (for a summary, see, e.g., Ref. 68 ). 
Testing these uncertainties, even the conservative 
choice of a 5 % uncertainty on the average mantle 
density is no problem for the Super-Kamiokande 
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detector. However, for the Hyper-Kamiokande 
detector, it should not be larger than about 2 %. 
Because of a partial averaging out of the uncer- 
tainties for the very long baselines considered in 
this paper, this uncertainty should be quite real- 
istic |5T] . 

Another source of uncertainty has to do with 
the knowledge of the location of a supernova, 
since it may affect the baseline length. In a future 
core-collapse supernova in our Galaxy, the elec- 
tromagnetic radiation can be obscured by dust 
in the interstellar space. In this case, the local- 
ization of the supernova can be done by study- 
ing neutrino-electron forward scattering. A con- 
servative estimate leads to restrict the supernova 
direction with an error of about 5° in the Super- 
Kamiokande detector 3.5, • However, with the in- 
crease of statistics expected in a megaton detec- 
tor, this uncertainty would be reduced. On the 
other hand, complementary information could be 
inferred by the triangulation method , namely 
the direct measurement of the delay of the sig- 
nal to the two detectors. Therefore, we consider 
that this small error would not imply a significant 
change in the long baselines that we are dealing 
with, since even for a 5° directional uncertainty 
the baseline length of about 12,700 km does not 
change very much. 

8. Summary and conclusions 

We have discussed the possibility to use the 
neutrinos from a future galactic supernova ex- 
plosion to obtain additional information on the 
Earth's core. First, we have summarized geophys- 
ical aspects and unknowns of the Earth's core. 
Then, we have investigated core-collapse super- 
novae as potential neutrino sources for a snapshot 
of the Earth's interior. Next, we have discussed 
the neutrino propagation from the production to 
the detection in detail, where we have especially 
focused on Earth matter effects on the neutrino 
oscillations of the supernova neutrinos. We have 
also put these effects into the context of the super- 
nova parameters, i.e., temperatures and deviation 
from energy equipartition. Furthermore, we have 
stressed the importance of supernova neutrinos 
arriving at the surface of the Earth as mass eigen- 



states for this technique, which we have also sup- 
ported by a discussion of decoherence of neutrino 
oscillations. We have shown that we need one 
detector on the surface of the Earth on the side 
towards the supernova, and another one in the 
shadow of the Earth's core. For the most likely 
scenario of not crossing the Earth's inner core, we 
have shown that the Earth's average core matter 
density could be determined up to 9 % with a 
Super-Kamiokande-like and 1.3 % with a Hyper- 
Kamiokande-like detector (all errors at the 2a 
confidence level). In addition, for a less likely two- 
parameter measurement of the outer and inner 
core matter densities, Hyper-Kamiokande could 
verify the existence of the inner core at the 3cr 
confidence level and measure the outer core mat- 
ter density with a precision of about 3.1 %, as well 
as the density jump between outer and inner core 
matter densities with a precision of about 59 %. 
The latter error is comparable to seismic wave 
geophysics, where, however, not the difference be- 
tween the average matter densities, but the mat- 
ter density jump at the outer-inner core boundary 
is measured. Thus, neutrino physics could pro- 
vide complementary information to geophysics. 
However, the quoted numbers for the precisions 
depend on the supernova parameter values, in- 
directly determined by the on-the-surface mea- 
surement, and are in some cases better, in others 
worse. Thus, the actual precisions will not be 
known before the supernova goes off. In general, 
we find that the more muon and tau neutrinos 
are produced, the larger the temperature differ- 
ence between the different flavors is, and the more 
different the degeneracy parameters for the differ- 
ent flavors are, the better the application works. 
Finally, we have discussed several uncertainties to 
these measurements and we have found that es- 
pecially the determination of the total muon an- 
tineutrino energy of the supernova causes prob- 
lems to our method. However, measuring not 
only electron antineutrinos, but also the other 
two flavors could completely eliminate the de- 
pendence on the supernova parameters. Further- 
more, the leading solar neutrino parameters have 
to be known with sufficient precision, which is 
about 0.2 % for Hyper-Kamiokande-like measure- 
ments. In summary, supernova neutrino tomog- 
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raphy could be a nice additional payoff of existing 
or planned detectors if all of the prerequisites can 
be met at the time when the next supernova ex- 
plodes. 
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